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Recently it has been proposed and experimentally demonstrated that a spin-orbit coupled multi- 
component gas in Id lattice can be viewed as spinless gas in a synthetic 2d lattice with a magnetic 
flux. In this letter we consider interaction effect of such a Fermi gas, and propose signatures in 
charge pumping experiment, which can be easily realized in this setting. Using 1/3 filling of the 
lowest 2d band as an example, in strongly interacting regime, we show that the charge pumping 
value gradually approaches a universal fractional value for large spin component and low filling of Id 
lattice, indicating a fractional quantum Hall type behavior; while the charge pumping value is zero if 
the Id lattice Filing is commensurate, indicating a Mott insulator behavior. The charge-density-wave 
order is also discussed. 


High spin quantum gas is a unique system of cold atom 
physics. The spin of atoms ranges from hyperfine F = 1 
of alkali atoms like Rb and Na, F = 9/2 ofK atoms 
to nuclear spin with SU(W) symmetry of alkali-earth- 
(like) atoms such as Yb and Sr, where W can be as 
large as ^ 10. Recently it emerges an interesting idea 
to use the internal spin degrees of freedom as another 
dimension, named as “synthetic dimension”, which nat¬ 
urally extends a D-dimensional system into a (D + 1 )- 
dimension one mm- In a Id lattice system, by applying 
two counter propagating Raman beams to couple differ¬ 
ent spin states, one can create a magnetic flux lattice 
in synthetic 2d [2]. This proposal requires a minimum 
amount of laser light and therefore minimizes heating 
from spontaneous emission. It also gives rise to a sharp 
edge in the synthetic dimension, which can help to visu¬ 
alize edge states. Very recently, two experimental groups 
have implemented this scheme, in Rb atom [3] and in 
Yb atom [4], respectively, and chiral edge states have 
been observed, for non-interacting (or weakly interact¬ 
ing) bosons [3] and fermions [4], respectively. Moreover, 
it is possible to create more exotic nontrivial geometry 
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The experimental setup and basic idea of synthetic di¬ 
mension are briefly illustrated in Fig. For instance, 
two Raman beams with tt and a polarization can couple 
spin state \m) to |m±l), where m can run from —F to F 
with totally W = 2F F 1 components. The Raman cou¬ 
pling has a spatial dependent phase factor where 

kji is the recoil momentum of Raman laser. We consider 
the situation 7 = 2/cRa = 27r/q (a is lattice spacing). The 
single particle Hamiltonian is therefore written as 

Ho = (1) 

j,m 

where j labels the site along the physical dimension x and 
m labels internal spin components, t is hopping along x 
and Cl is Raman coupling strength. There are two differ¬ 
ent but equivalent pictures of this single-particle Hamil¬ 
tonian 


ofWW\o 

IpuS 

FIG. 1: Illustration of two different pictures of this system 
and the idea of charge pumping, (a) Two Raman beams along 
X are applied to a multi-component quantum gases in opti¬ 
cal lattices. The Raman beams couple different spin states 
and generate a coupling between spin and momentum kx ■ (b) 
Different spin components are viewed as another dimension, 
the system is mapped into a 2d spinless particle in a mag¬ 
netic held. By applying an electric held along the physical 
dimension x, it generates a charge pumping along the syn¬ 
thetic dimension, which can be detected by measuring spin 
populations. 

(a) A Id system of high spin atoms with spin-orbit 
coupling: by applying a spin and site dependent rotation 

^ the Hamiltonian Eq. [^becomes 

Ho = + h.C.) (2) 

j,m 

The spin dependent hopping term, together with a con¬ 
stant spin flipping coupling, gives rise to the spin-orbit 
coupling effect, which have been extensively discussed for 
spin-1/2 case in continuum in the past a few years [6]-[9]. 

(b) A 2d system with magnetic flux: If we view m as 
another synthetic dimension, say, labelled by y, Eq. 
represents a situation that spinless atoms hops in a 2 d 
space, with hnite number (VF) of chains and open bound¬ 
ary condition along y. More importantly, hopping along 
a close loop around a plaquette accumulates a phase of 
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e*'^, that is equivalent to say, each plaquette has a flux of 
(^0 is magnetic flux unit), which is gauge coupled 
to motion of atoms. 

These two equivalent pictures build up an intriguing 
connection between spin-orbit coupled high spin particles 
in Id and spinless but charged particles in a 2d ladder 
geometry with magnetic flux. For instance, at single¬ 
particle level, the chiral edge current from the picture 
(b), as observed in Ref. [3j |4], is equivalent to spin- 
momentum locking effect from picture (a), as observed 
in spin-1/2 case [^-[8]. In this letter we aim at study¬ 
ing effect of repulsive interaction in this system [10] . For 
2d electron gas in a magnetic field, it is known that frac¬ 
tional quantum Hall (FQH) state will emerge with strong 
repulsive interactions. However, in this system there are 
a few important differences that worth emphasizing first: 

1. In the synthetic dimension, the system always only 
has finite number of chains. Though it is possible to cre¬ 
ate a periodic boundary condition in synthetic dimension 
with some complicated laser setting, in the most simple 
and natural setup, the system has an open boundary con¬ 
dition in synthetic dimension. Thus, the finite size effect 
can be significant. 

2. Along the physical dimension, the interaction is on¬ 
site and short-ranged. While in the synthetic dimension, 
the interaction is long-ranged. Let us consider a SU{W) 
invariant interaction m, 

-^int — U ^ ^ '^j,m'^j,m' 5 (3) 

atoms in any two sites along the synthetic dimension, 
despite of their separation, interact with the same inter¬ 
action strength. In another word, in the 2d lattice the 
interaction is very anisotropic. 

3. For normal FQH effect, the only relevant parameter 
u is the ratio between fermion number to flux number. 
In our case, it will be 

= J^= ^ ^ (A) 

Afflux WL/q WL^ ^ ’ 

where N is total number of fermions, L is the number 
of sites along x-direction. However, from the picture (a) 
that our system is one-dimensional, it is natural to intro¬ 
duce another filling factor 


i^id = ^, (5) 

and if is an integer, one may expect a trivial Mott 
insulator rather than a FQH state when interaction U is 
sufficiently large. Thus, when u is fixed, we still have two 
other tunable parameters, i.e. z^id and W. 

Hereafter we shall fix = 1/3 as a typical example. 
Given those differences mentioned above, one may won¬ 
der whether one will still have a FQH-type behavior un¬ 
der strong repulsive interactions. Beside, whether there 



FIG. 2: Charge pumping for different z^id and W. Numbers 
in box are charge pumping Q after insertion of one flux. Blue 
triangles are points calculated by ED and red circles are points 
calculated by DMRG, with u fixed at 1/3, Q = t and U = 6t. 
(a-c) marked three cases where spectral flow will be shown in 

Fig. 13 


is an effective scheme to detect such a state in this cold 
atom setting. The rest of this paper is devoted to answer 
this question. 

Charge Pumping. To reveal the interaction effect, we 
propose to perform a charge pumping experiment utiliz¬ 
ing the advantage of synthetic dimension. Let us consider 
applying an electric field along x direction. In our numer¬ 
ical calculation below, this is implemented by a periodic 
boundary condition in x direction, i.e. T(I/) = e*^T(0). 
In cold atom experiment, this can be realized by, for in¬ 
stance, applying a field gradient or moving lattice with 
a constant velocity. Gharge pumping here means charge 
transfer along the synthetic dimension [12], as schema¬ 
tized in Fig. 13 Defining the charge center as 

j,m 

the charge pumping after inserting one flux is given by 

Q = Y{0 = 27t)-Y{0 = 0). (7) 

If it is in real space, detecting charge transfer requires in- 
situ image, while charge transfer in synthetic dimension 
means changing of spin population, which can be much 
easily detected from the Stern-Gerlach experiment in the 
time-of-flight image [nun]. 

In Fig. 13 we present the charge pumping value Q 
for various z/id and VF, with u fixed at 1/3. This result 
is obtained by numerically solving the many-body wave- 
functions with Hamiltonian H = Hq + ^int? either by 
exact diagnolization (ED) or density matrix renormal¬ 
ization group (DMRG) methods. For ED the maximum 
number of particle is six and the dimension of Hilbert 
space is of the order of 3 x 10^. Eor DMRG, the maxi¬ 
mum number of particle is ten, and the truncation error 
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FIG. 3: Left column: Spectral flow under the insertion of flux 
(i.e. changing periodic boundary condition 0 from zero to 
Gtt.) Right column: charge center y as a function of 0. (a-c) 
correspond to different W and uid as marked in Fig. (a) 
W = 9, uid = I {N = L = 5); (b) VF = 5, uid = 5/6 {N = 5, 
L = 6) and (c) W = 14, = 2/3 (N = 4, L = 6). All these 

cases are calculated by ED. 


is of the order of 10“^. Each eigenstate has a well-defined 
quantum number K that is the center-of-mass momen¬ 
tum along X. We also plot how these eigenstates evolve 
under the changing of 6>, as shown in Fig. 3, and for 
ground state, we calculate F as a function of 0 with Eq. 
[^and deduce Q with Eq. 

We find following features: (i) For uid = 1, Q is identi¬ 
cally zero for all W. In this case there is always a unique 
ground state which will not interchange with other states 
under flux insertion, as shown in Fig. 3(a). This Mott 
insulator phase with commensurate uid is due to one¬ 
dimensional nature of the system, or equivalently to say, 
anisotropic nature of interaction in picture (b). (ii) Since 
here we consider q is an integer, the smallest value for q is 
q = 2. With q = 2 and = 1/3, uid = vWjq can at most 
be LF/6. We note that for g' = 2, the Hofstadter spec¬ 
trum exhibits Dirac cone instead of fully gapped band, 
and the lowest band does not have well defined Chern- 
number. Therefore, for v\d = WI^ (with LF ^ 6) we also 
find Q = 0. A typical spectral flow and charge pumping 
is shown in Fig. 3(b). (hi) For small W, or for larger FF 
but vid closer to unity, though Q is generally non-zero, 
it takes non-universal values. This fluctuating Q indi¬ 
cates some Fermi-liquid type states m- (iv) For large 
FF and smaller z^id, Q gradually approaches a universal 
fractional value of 1/3. A typical spectral flow under flux 
insertion is shown in Fig. 3(c). One can see that there 
are three low-lying states (though not exactly degener¬ 
ate) exchange one with the other as 0 increases, and the 



FIG. 4: Gharge pumping Q as a function of Raman coupling 
strength Q/t. The dashed line is 1/3. 


spectrum recovers itself only after 0 changes Gtt. These 
features are consistent with a FQH effect. This is because 
for large FF, the finite size effect in synthetic dimension 
becomes insignificant, and for smaller uid away from com¬ 
mensurate filling, the lattice effect also becomes weaker. 
Previously, it has been shown by Luttinger liquid theory 
that for continuum models, fractional state can emerge in 
a one-dimensional system with spin-orbit coupling m- 

We also study the charge pumping value Q as a func¬ 
tion of Raman coupling strength D, as shown in Fig. 
1^ As the synthetic magnetic field is resulted from Ra¬ 
man coupling, one naturally expects that Q will vanish 
as D ^ 0. Indeed, we show in Fig. |^that when 0/t is 
smaller than certain value, Q starts to derivates from 1/3 
and fast drops to zero. This feature is particularly clear 
for large FF (e.g. green points for FF = 14 in Fig. [^. 
We have also looked at Q for smaller U/t and find when 
U < t, Q also takes fluctuating non-universal values. 

Periodic Boundary Condition in Synthetic Dimension. 
We also find that if one applies a more involved laser 
setting to achieve a periodic boundary condition along 
the synthetic dimension, it will greatly help to stabilize a 
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FIG. 5: Energy diagram with periodic boundary condition in 
synthetic dimension, (a) Energy levels for different momen¬ 
tum Ky -h WKx. (b) Energy of the lowest three states as a 
function of 0. Here VF = 4, uid = 1/3 {N = 4, L = 12). This 
case is also calculated by ED. 
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FIG. 6: Fourier transformation of density p{qx) for different 
W and z^id 

fractional state. For instance, for the cases with VF = 4 
we presented in Fig. we do not find accurate fractional 
charge pumping. While when we apply a periodic bound¬ 
ary condition in synthetic dimension, we show the energy 
level for different momentum (labelled by Ky + WKx as 
now both Kx and Ky are good quantum numbers) in Fig. 

We find a very accurate three-fold degeneracy, with 
energy splitting smaller than and these states are 

separated from other excited states by a gap ^ 0.04t. The 
total momentum of these three states are also consistent 
with generalized Pauli-exclusion principle analysis m- 
Moreover, these states exchange one and the other under 
the flux insertion, and do not intersect with other excited 
states, as shown in Fig. ib). By calculating Berry cur¬ 
vature with twisted boundary conditions in both physical 
and synthetic dimensions m, we numerically find that 
their many-body Chern number Ci = C 2 = Cs = 0.333. 

Density Order. Finally we look at real-space charge 
density-order. We consider the onsite total density pj = 
M we plot the Fourier transform of pi 
as p{qx) = (1/T) (p = is the average 

density). p{qx) shows a clear peak at qx 1 ( 211 ) = and 
qxl(2i\) = 1 — z^id- This feature exists for both open 
and periodic boundary condition in synthetic dimension. 
Similar situation has also been found in several other 
models m- This is reminiscent of usual FQH state in 
the thin torus limit m- 

Conelusion. In summary, we have studied interaction 
effects in the synthetic dimension picture of high spin lat¬ 
tice Fermi gases with Raman-coupling induced spin-orbit 
coupling. Our studies are mainly focused on signatures of 
charge pumping experiment, which becomes much easier 
in this setting, as charging pumping along the synthetic 
dimension can be visualized by measuring spin popula¬ 
tion. For fixed = 1/3 case, we investigated how the 
charge pumping value depends on number of spin com¬ 
ponent VF, fermion density z^id, and Raman coupling fl/t. 
We conclude that a universal fractional charge pumping 
Q = 1/3 is favorable for VF ^ 1, z/id ^ 1, 0/t ^ 1 and 
with strong interactions. Similar results have also been 


obtained for strongly interacting bosons. 
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Note Added. When finishing up the manuscript, we 
become aware of another paper Ref. [20] in which the 
same system is studied by DMRG. Charge density wave 
order is also discussed in this paper. 
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